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ABSTRACT 
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New visible polarization data combined with existing IR and FIR polarization 
00 ■ data are used to study how the magnetic field threading the filamentary molecu- 

lar cloud GF 9 connects to larger structures in its general environment. We find 
that when both visible and NIR polarization data are plotted as a function of 
extinction, there is no evidence for a plateau or a saturation effect in the polariza- 
tion at Ay « 1.3 as seen in dark clouds in Taurus. This lack of saturation effect 
suggests that even in the denser parts of GF 9 we are still probing the magnetic 



field. The visible polarization is smooth and has a well-defined orientation. The 
IR data are also well defined but with a different direction, and the FIR data in 



the core region are well defined and with yet another direction, but are randomly 
distributed in the filament region. On the scale of a few times the mean radial 
dimension of the molecular cloud, it is as if the magnetic field were 'blind' to the 
spatial distribution of the filaments while on smaller scales within the cloud, in 
the core region near the IRAS point source PSC 20503+6006, polarimetry shows 
a rotation of the magnetic field lines in these denser phases. Hence, in spite of 
the fact that the spatial resolution is not the same in the visible/NIR and in the 
FIR data, all the data put together indicate that the field direction changes with 
the spatial scale. Finally, the Chandrasekhar and Fermi method is used to evalu- 
ate the magnetic field strength, indicating that the core region is approximately 
magnetically critical. A global interpretation of the results is that in the core 
region an original poloidal field could have been twisted by a rotating elongated 
(core+envelope) structure. There is no evidence for turbulence and ambipolar 
diffusion does not seem to be effective at the present time. 
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1. Introduction 

With the measurement of visible linear polarization of field stars, followed by develop- 
ments of Near and Far Infrared (respectively NIR and FIR) and submillimetric polarimetry, 
it is now possible to trace magnetic field structures from the scale of the Milky Way to that 
of star forming regions and cores. In the diffuse Interstellar Medium (ISM), the general 
assumption that elongated dust grains usually have their smallest axis aligned in a direction 
comparable to that of the local magnetic field has lead to a toroidal representation of the 
Galactic magnetic field. A similar topology of the magnetic field can also sometimes be 
observed in other galaxies (e.g. Scarrott et al. (1990)). On the smaller scales of molecular 
clouds, a similar relation exists between grain alignment and magnetic field directions (see 
Lazarian (2003) for an exciting review on the subject) and many topologies of the mag- 
netic field are suspected to exist. Among others, helicoidal, toroidal, poloidal and hourglass 
configurations have been proposed to interpret observed polarization patterns, taking care 
of the physical conditions in these environments (e.g., Orion A, M17). However, one of the 
most important limitations of linear polarimetry is that we only get a two-dimensional pro- 
jected average representation of the magnetic field topology because of the integration on 
the plane of the sky of the contribution to the polarized emission due to each aligned grain, 
or its dichroic absorption for visible/NIR observations. This leads to a degeneracy problem 
since many three dimensional (3D) topologies of the field can produce identical polarization 
patterns. A method to get rid of this degeneracy problem is the use of models. On the scale 
of relatively nearby filamentary molecular clouds and cores, this method has been used and 
allows in some cases to select one field topology among many possible ones (e.g., Matthews 
et al. (2001), Vallee et al. (2003)). 

In addition to models, other diagnostics of magnetic fields can be used to lift the degener- 
acy due to the plane-of-sky integration. When column densities are high enough, observations 
of the Zeeman effect can yield the line-of-sight component of the magnetic field strength (e.g. 
Crutcher (1999)). Finally, by determining the ion-to- neutral molecular line- width ratio in 
non turbulent regions with a linear flow, it is possible to determine the angle between the 
field direction and the line of sight (Houde et al. (2000a), Houde et al. (2000b), Houde et 
al. (2001)). By combining all three methods together, we get the 3D field topology (e.g., 
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Houde et al. (2002) for M17, Houde et al. (2004) for Orion A). The observation of Zeeman 
splitting of spectral lines requires bright dense cores, but even in other regions, it is possible 
to get the 3D geometry of the field even if the field strength is unknown. 

These methods can inform us about how magnetic fields can contribute to cloud dynam- 
ics, but our ultimate goal is to understand the competition between gravitational, thermal, 
magnetic, and turbulent forces. More specifically, magnetic fields and turbulence are both 
competing to slow down the star formation process, otherwise stars would form much faster 
than observations show. Combining polarization modeling techniques (e.g., Fiege & Pudritz 
(2000)) and predictions for magnetized simulations of turbulence (e.g., Heitsch et al. (2001), 
Padoan et al. (2001)), one can hope to be able to determine the relative contributions of 
magnetic fields and turbulence. 

To reach this goal, and to understand how magnetic fields thread the dark filaments, and 
to put observational constraints on modeling we propose to make comparisons of magnetic 
fields as seen in and around individual objects. This approach, namely a multi-scale analysis 
of the magnetic field, has begun to be done via comparisons of visible, NIR, FIR and/or 
submm polarization data in the OMC1 (Schleuning (1998)), OMC3 (Matthews et al. (2001)) 
and GF 9 regions (Jones (2003)). In these works, aspects of the galactic magnetic field 
extending on tens of degrees are compared with aspects of the clouds probed in the densest 
regions of the clouds with resolutions ~ 10", but there is generally a lack of information 
between these two scales. Here we present observations at an intermediate scale, as much 
as the presence of field stars allows it, to have a clearer view of the topology of the field at 
the interface ISM/filamentary molecular clouds. From a general point of view, this approach 
should provide new insights on the way magnetic fields and clouds are dynamically coupled 
and shaped. In this paper we address this problem for the GF 9 filamentary molecular cloud 
region. A forthcoming paper will address the more complex situation for the Orion A region. 

A description of data available about GF 9 is given in section 2. Observations and 
observational results are presented in section 3. An analysis and a discussion follow in 
section 4, and our conclusions are given in section 5. 

2. The GF 9 Region 

The GF 9 dark filamentary structures have been catalogued first as B150 by Barnard 
(1927), L1082 by Lynds (1962) and GF 9 by Schneider & Elmegreen (1979). A first estimate 
of the distance to these condensations as d = 440 pc was determined by Viotti (1969) (See 
Hilton & Lahulla (1995)) but an estimate of d = 200 pc to the extreme class object GF 
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9-2 based on star counts was reported more recently by Wiesemeyer (1998). The first 
estimate is reported in catalogs or works on relatively large core samples by Benson & Myers 
(1989), Goodman et al. (1993), Dobashi et al. (1994), Lee & Myers (1999), Furuya et 
al. (2003), while the second one is reported in works by Andre et al. (2000), Furuya et al. 
(2003) and Froebrich (2005). Assuming d = 200 pc, GF 9-2 (see Figure 1) has a luminosity 
-^boi = 0.3L Q , an envelope mass M env ~ O.5M , a ratio £ S ubmm/£boi ~ 10% an d a bolometric 
temperature T bol < 20K. No outflow manifestation, nor structures are detected but infall 
motions are observed (Glisten (1994), Wiesemeyer (1997), Wiesemeyer (1998) and Andre 
et al. (2000)). GF 9-2 must be particularly young since it has a H 2 maser but no detected 
outflow (Furuya et al. (2003), Furuya et al. (2005)). 

Two ?y 8' x 10' regions, namely GF 9-core and GF 9-fila, have been observed by Ciardi 
et al. (1998), Ciardi et al. (2000). The two areas are shown in Figure 1 by blue boxes 
identified by C and F. The GF 9-core region is associated with the IRAS point source class 
protostar PSC 20503+6006, but also with GF 9-2 and an object named [LM99]351. JHK 
photometry of these two regions probe similar masses but the filament region is consistent 
with a uniform-density cylindrical cloud, while the core region contains a centrally condensed 
core with the highest extinction region lying to the north of the IRAS point source. CO, 13 CO 
and CS surveys of these two regions give consistent views with these two pictures, but also 
reveal infall/outflow motions in order to explain CS line broadening in the core region where 
n H2 « 15000 ± 3000 cm" 3 (assuming d = 440 pc). In the filament region, n H2 « 6000 ± 1200 
cm -3 (see Table 2 in Ciardi et al. (2000)). CO and NH 3 observations of three regions with 
an infrared association are presented by Benson & Myers (1989). These sources are IRAS 
20520+6003, 20526+5958 and 20503+6006, and they are associated respectively with the 
LDN 1082A, B and C cores. An analysis of these data by Goodman et al. (1993) shows 
that motions consistent with uniform rotation in dense cores are present in LDN 1082A and 
LDN 1082B, while not in LDN 1082C. However recent works by Furuya et al. (2005) show a 
velocity gradient along the major axis in the core in GF 9-2 which could be due to a rotation 
of the core. 13 CO observations by Dobashi et al. (1994) have also been done at 3 positions 
in GF 9. Referring to Lee & Myers (1999), objects named [LM99]349 and [LM99]350 (see 
figure 1) have no associated embedded young stellar objects (EYSO) while [LM99]351 1' 
west of GF 9-2) and [LM99]352 have associated EYSO. 

First A = 184 /xm FIR polarimetry observations from space with ISO in GF 9 were 
done in the core and filament regions by Clemens et al. (1999). Additional measurements 
at A = 1.65 /im on 9 stars shining through the same regions of GF 9 are also reported by 
Jones (2003). Finally, 850 /iin polarimetry observations were done at the JCMT by Jane S. 
Greaves in 1998 at one position in the filament. 
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Initially, the relatively small heliocentric distance of GF 9, its filamentary shape, and 
potentially interesting FIR measurements in this region were the major factors which moti- 
vated us to conduct visible polarimetry around these filamentary dark structures. 

3. Observational results 
3.1. Observations 

The observations were carried out on the 1.6 m telescope at the Observatoire du Mont- 
Megantic (OMM), Quebec, Canada, between 2000 September and 2003 July using a 8.2 
arcsecond aperture hole and a broad band red filter (RG645: 7660 A central wavelength, 
2410 A FWHM). Polarization data were taken with Beauty and the Beast, a two-channel 
photoelectric polarimeter, which uses a Wollaston prism, a Pockels cell, and an additional 
quarter wave plate. The data were calibrated for polarization instrumental efficiency, instru- 
mental polarization (due to the telescope mirrors) and zero point of position angles using 
a prism, non polarized and polarized standard stars, respectively. On average, the instru- 
mental polarization was 0.081 ± 0.061% and was subtracted. The observational errors were 
calculated from photon statistics and also include uncertainties introduced by the previously 
mentioned calibrations. The final uncertainty on individual measurements of the polariza- 
tion P is usually around 0.1%. For more details on the instrument and the observational 
method, see Manset & Bastien (1995). 

The JCMT data archive has been used and a treatment of the data taken in 1998 was 
done with commonly used SCUBA data reduction routines. Data were taken in relatively 
strong water vapour absorption (rcso > 0.08) conditions, and we finally got very low signal 
to noise polarization ratios and no significant results out of noise were found. 

3.2. Visible polarization measurements 

Most of the stars were chosen with the GSC interface at the Canadian Astronomy Data 
Center (CADC). The others were selected directly on Digitized Sky Survey (DSS) red plates. 
All stars observed are compiled in Table 1. Column 1 gives the GSC number. Column 2 
shows a GF 9 star 1 designation number used in this paper. Columns 3 and 4 give equatorial 



1 We use an S before star numbers in GF 9 to distinguish them from other numbers already in use to 
designate sources detected at longer wavelengths. 
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coordinates at epoch 2000. Degrees of polarization and equatorial position angles with their 
uncertainties are given in columns 5 and 6 respectively. Column 7 gives the visible extinction 
coefficient Ay from the Dobashi et al. (2005) atlas. Finally, 1 is tabulated in column 8 
when P > 3op (meaning o$ < 9.5°), otherwise is tabulated in this column. We note that 
when a star has no GSC number, not much information is generally available about it. 

Figure 2 presents a polarization map of stars for which the polarization degree meets the 
condition P > 3<7p. This map is to be compared to Figure 1. The histogram of the degree of 
polarization is shown in Figure 3. If we only consider data for which P > 3<jp (data shown 
with full lines), the mean and the standard deviation of the distribution are P = 2.22% and 
Sp = 1.33%, respectively. The histogram of the position angles is shown in Figure 4. If we 
only consider data for which oq < 9.5° (data shown with full lines), the mean position angle 
and the standard deviation of this well peaked distribution are 9 = 127.4° and S$ = 25.5°, 
respectively. This means a standard deviation of 0.445 rad, comparable to those found for 
clouds with embedded clusters (see Myers & Goodman 1991b). 

Stars S35, S70, S75, S86 and S92 are identified on both maps in Figures 1 and 2. Stars 
S35 and S86 are suspected to have variable linear polarization. These stars are located in the 
continuity of two filaments but in gaps where extinction appears to be relatively low. The 
data relative to these stars are compiled in Table 2 with the respective Julian Date (J.D.) of 
each measurement appearing in the last column. 

The position angle of stars S70, S75 and S92 differ significantly from the mean direction 
of the other stars (9 = 127 A with So = 25.5°). Star S70 is a visual binary whose components 
fit both in the aperture hole used for the observations; its position angle is perpendicular 
to this mean position angle. Thus, we may be in the presence of a physical binary system 
with circumstellar material, or of a single system presenting intrinsic polarization and a 
background star. The difference of the position angles of stars S75 and S92 with the average 
value may be related to their position relative to the clouds as seen on Figure 1. Star S75 
(Ay = 3.3) is located north-east of IRAS 20520-6003 at the edge of two filamentary structures 
mainly parallel. Comparatively, star S92 is also located in a dense region (A v = 2.1) in a 
small filamentary structure oriented mostly north-south. Given their large extinction, these 
two sources could either be deeply embedded in their clouds or be young stellar objects 
surrounded with disks in which multiple scattering can produce a relatively high degree of 
polarization (e.g., Bastien & Menard (1990)). Without multiwavelength data it is difficult 
to distinguish between these possibilities. 

Three stars in our sample, GF 9 SI, GF 9 S2 and GF 9 S15 have measured spectral 
types and visual magnitudes. However we could not get reliable distances from these values. 
The Ay values obtained are negative or yield very large distances. Similarly using Ay values 
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from Dobashi et al. (2005), we could not get reliable distances. 

3.3. Variations of P with A v 

In order to study the variations of the degree of polarization with visible extinction, 
we used the Dobashi et al. (2005) catalog to determine Ay, in the direction of each star 
for which polarization measurements are available. This catalog is the first version of the 
atlas and catalog of dark clouds derived by using the optical database Digitized Sky Survey I 
(DSS), and applying a traditional star-count technique to 1043 plates contained in the DSS. 
Using the visible extinction coefficient Ay listed in Table 1, values of P(A = 7660A) against 
E-q-y, where E-q, — V — 3~l"' 

shown in Figure 5. For all stars, except GF 9 S35 and GF 9 
S86 (suspected variables), data directly observed at this wavelength are shown with crosses. 
Assuming NIR data from Jones (2003) and visible data probe the same types of grains, the 
Serkowski law, 

P(A) = P ma ,exp[-Kln 2 (^p)], (1) 

where K = 1.15, is used to derive the relation Py = 1.97Ph- This relation is used to 
convert the H-band data to equivalent A = 7660A data, assuming A max = 7660A. These 
equivalent data are listed in Table 3 and are shown with diamonds in Figure 5. As was 
previously mentionned by Jones (2003), no saturation of P can be seen at about Ay « 1.3 
(or Eb-y ~ 0.4), and Py ~ 1.6%, as observed in cold dark clouds in Taurus (See Arce et al. 
(1998)). Data in the visible confirm this analysis but over the whole filament and the region 
surrounding it, and not just in the core region. Thus, while not universal, it is possible in 
some regions to probe magnetic fields inside dark clouds via polarimetry of background stars 
as suggested here by data taken in the visible and in the NIR. 

4. Analysis and Discussion 

4.1. Ambient Magnetic Field Orientation in the Vicinity of GF 9 

As can be seen in Figure 2 and in the histogram of position angles in Figure 4, the 
polarization pattern in the vicinity of GF 9 is relatively well oriented on the plane of the sky 
with a mean direction of 9 = 127 A and a standard deviation So = 25.5°. Since the filament 
has an incurved concentric shape reminiscent of supernova remnants (Schneider & Elmegreen 
(1979)), one can ask if a possible object at the origin of this shape could be located somewhere 
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in the region delimited by the black circle traced in Figure 1. A possible candidate seems 
to be the post AGB star GLMP 1012 (see Figure 1) but given its distance estimate d = 3 
kpc (Preite-Martinez (1988)) compared with those of GF 9, even by considering 50 km s _1 
winds, it is difficult to defend the view that this source can be responsible for the shape of 
the filaments. However since the problem is tridimensional, a source at the origin of such 
a shape could also be located outside of the domain defined by this circle. Maybe such a 
source exists but has not been detected, or the apparent concentric shape of the dense ISM 
is only a projection effect and other mechanisms, such as turbulent fluctuations in the ISM 
or density waves, are at the origin of its formation. A quick look in the catalogs reveals no 
potential candidate in a circle a few degrees in diameter around the position R. A. (2000) = 
20 h 50 mn 53 s and Dec. (2000) = 59°51'59". 

If now one compares the direction of the inferred magnetic field with the orientation of 
the filamentary structures on the plane of the sky, it seems that the magnetic field has 
approximately the same orientation over an arc circle which covers about 135°, from a 
position angle (P.A.) = 270° to 45° as measured from the center of the circle in Figure 
1. When looking at both maps from east to west in Figure 1 and 2, a slow rotation of the 
direction of the field is present, the average directions near the core and near the filament 
regions differ by about 10° (see Table 4), significantly under the variation in P.A. of the 
filament over the whole region. Moreover from east to west, there is a great dispersion 
in the orientation of the filamentary structures, thus it is hard with this two dimensional 
representation to see a direct impact of the field on the apparent shape of the filaments. This 
is in agreement with the fact that in other regions there is no prefered projected angle on 
the plane of the sky between the directions of the ambient magnetic field and the orientation 
of the filamentary structures (e.g. Myers & Goodman (1991b), Heyer et al. (1987)). See 
also Heiles et al. (2000) for a discussion on this subject. 



4.2. Multiscale Analysis of the Magnetic Field 

4-2.1. From filamentary clouds scale to cores scale 

Comparisons of the orientation of the magnetic field according to the visible, IR and 
FIR data around and in the core and the filament regions can be done by looking at the 
maps in Figures 6 and 7 respectively. In the core region, when going from visible to IR, to 
FIR data, one can see an anticlockwise rotation of the P.A. by about 63° in total. For each of 
these two regions, the means and dispersions of the degree of polarization and of the P.A. of 
the visible, IR and FIR data are compiled in Table 4. The P.A. for the FIR data corresponds 
to the observed polarization, rotated by 90° to give the direction of the magnetic field. In 
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Figures 6 and 7, only FIR data for which the P.A. is approximately the same when using the 
KLC instrumental polarization and the C-off instrumental polarization (Clemens 1999) are 
shown. The degree of polarization of the FIR data is not considered in the present analysis. 
Table 4 also includes similar information for the whole GF 9 region and for the local galactic 
scale, to be discussed below. 

Comparison of Visible and IR data 

First, visible and IR data do not probe the same regions on the plane of the sky. Visible 
data probe the more diffuse parts around prestellar cores while the IR data probe the denser 
phases where protostellar cores are present. Secondly we note that in the core region the IR 
vectors have approximately the same orientation which implies that the IR polarization is not 
intrinsic to each star but produced by dichroic absorption by aligned grains. This statistical 
argument, however, does not apply in the filament region where only one IR measurement 
with sufficient signal-to-noise ratio is available. Thirdly, in both cases the resolution used 
with each technique is comparable. Fourthly, the simple fact that the 'IR' magnetic field 
has not the same orientation than the 'visible' one implies that if the grains observed with 
these two techniques have identical properties, and if their dichroic absorption properties are 
the same at both IR and visible wavelengths, then the net polarization produced by aligned 
grains located in the densest parts of the cloud dominates the net polarization produced 
by aligned grains located in the foreground diffuse ISM. All these facts, in addition to the 
dependence of P with Ay mentionned in section 3.3, suggest that there is a real rotation of 
the magnetic field lines in the core region when moving across the plane of the sky. 

Comparison of IR and FIR data 

Firstly, IR and FIR data probe the same regions on the plane of the sky in both the core 
and the filament regions. Secondly we note that in the core region the FIR vectors are very 
well aligned but not with the same P.A. than the IR vectors. IR data result from a magnetic 
field oriented at P.A. ps 171° while FIR data suggest a magnetic field oriented at P.A. ps 19°. 
Thirdly, the resolution with these two techniques is not the same. The FIR data resolution 
is 80", many times the 'pencil' beam resolution of the IR and visible data. Finally, FIR 
observations probe the radiation emitted by grains mainly located in the densest parts of 
the cloud while IR observations probe grains located between the stars and the observer. 

When looking at Figure 6, we see that IR and FIR data cover a common area ~ 4' x 4' in 
size. Thus if these two observational techniques probe the same column densities but at two 
different resolutions, one should be able to explain the rotation of the magnetic field lines as 
seen at these two wavelengths. A first hypothesis is that a given magnetic field morphology 
leads to several orientations according to the resolution at which it is observed (See Fiege 
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& Pudritz (2000), but also Heitsch et al. (2001)). However, here we see that IR data 
are spatially distributed through FIR data, thus by using a simple principle of juxtaposition 
it seems impossible to reproduce the mean orientation of FIR vectors seen through m 80" 
beams by combining several IR vectors with very small cross sections. On the other hand, 
since the mean position angle of the stars observed in the visible is not the same than the 
one for stars observed in the IR, these IR-observed stars should not be foreground to the 
clouds, but rather embedded in or background to the clouds. In what follows we review 
possible explanations for the shift in position angles. 

(1) The stars observed in the IR could be located inside the GF 9 core. If this was 
the case, the impact of their presence on their immediate environment should be observed. 
IRAS observations of these stars show that their luminosities satisfy L suhrnm / L ho \ ~ 10% 
(see section 2), which would argue against them being Class sources, i.e., very young. 
Therefore, we reject the hypothesis that these IR stars are all embedded in the GF 9 core. 

(2) The stars observed in the IR could be background to the GF 9 core region. This is 
suggested by their position in the P — £b-v diagram shown in Figure 5. Two possibilities 
can be invoked to explain the difference in position angles: 

(a) A first and may be naive approach is assume that grains along the line of sight to the 
GF 9 core are all at the same temperature, and that both the IR and FIR observational 
techniques probe all grains with the same efficiency If this is right, the mean position angles 
observed at both wavelengths should be similar, except for the known 90° difference between 
them. In this scenario, the difference in position angles is explained by a second cloud or 
filamentary structure located behind the GF 9 core, and in which dust grains would be 
aligned differently. 

(b) A second and probably more realistic scenario is to emphasize differences related to the 
intrinsic nature of both observational techniques. Based on point (1) mentioned above, the 
assumption of a dust temperature gradient for an externally heated cloud seems reasonable. 
In this case, the FIR measurements should be more sensitive to "cold" dust grains aligned 
in the densest regions of the core than to "hot" dust grains aligned in the envelope and 
in the diffuse ISM. On the other hand, IR observations should probe all the grains along 
the line of sight, independently of their temperature. Thus grains aligned differently in the 
envelope than in the denser core should both contribute to the IR polarization position angle, 
explaining the offset in position angles. 

While scenario (2b) does not reject the possibilty of a second cloud or filament evoked 
in (2a), it has the benefit to propose a simple and consistent explanation. Thus we prefer 
scenario (2b), in which the intermediate position angles of IR vectors shown in figure 6 may 
result from a vectorial addition of both FIR and visible vectors covering this region and its 
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vicinity. 

Comparisons of the magnetic fields orientations with the cores elongations 

We compare in Table 5 the orientations of five protostellar cores spread into the fil- 
aments with the magnetic field orientations. Orientations of the cores L1082A, B and C 
were determined with the (J, K) = (1,1) lines of NH 3 while orientations of the LM cores are 
based on the visible DSS data 2 At first sight we see that, except for L1082B, these dark 
condensations have a tendency to be elongated along the filaments as was noticed by Myers 
et al. (1991a) (see Figures 1 and 2). Names, positions and P. A. of these cores are shown in 
Table 5 in columns 1 through 4 respectively. For objects LM351 and L1082C, located in the 
core region, and for objects LM349 and LM350, located in the filament region, the offset, 
|A$7660 1, on the plane of the sky between the P.A. of each protostellar core and the mean 
magnetic field orientation probed with visible data is shown in column 5. The same is done 
with the IR and FIR measurements in columns 6 and 7 respectively. These offsets were not 
estimated for L1082B since the magnetic field direction is undersampled around this object. 

On larger scales, the magnetic field is mostly along the minor axis of the cores, as can 
be seen in Figure 2. However, on smaller scales, if we accept that IR measurements can 
probe the magnetic field in dense regions as it is suggested by the variations of P with Ay 
shown in Figure 5, we see that a rotation of the field is apparent in the core region. On 
this smaller spatial scale the field is well aligned in the core region in a direction mainly 
parallel to the cores' major axes. Turbulence models such as those developed by Padoan 
et al. (2001) show that matter may flow onto cores along filaments with the field being 
stretched along the filaments. In such a case, the field would tend to be along the major 
axis of the elongated low-density structures connected with higher density cores. However, 
such models show that the polarization pattern probing dust in the low density structures 
would be chaotic rather than smooth and well aligned. Moreover, Furuya et al. (2005) 
showed recently that a velocity gradient of 2.3 km s -1 pc _1 is present along the major axis 
of the core GF9-2 but none along its minor axis. Therefore, if we assume that sources named 
LM351, L1082C and GF 9-2 are associated to a common spatial volume traced with different 
tracers, and if we assume that the velocity gradient along the major axis of GF 9-2 is due to 
rotation of the core, this picture suggests that an original poloidal field has been twisted on 
small scales by the rotation of the core into a toroidal field. This would be consistent with 
a magnetic support model. 



2 In fact, objects L1082C and LM351 may be the same source. The difference in coordinates and in P. A. 
comes about because of the different observational techniques used, NH 3 for L1082C and visible extinction 
for LM351. 
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On the other hand, this effect does not seem to be present in the filament region where 
the magnetic field probed with IR data is mainly parallel to the small axes of the cores 
detected in this area, in approximately the same direction than the magnetic field seen at 
greater scales. But on the other hand if FIR measurements are reliable, the chaotic FIR 
polarization pattern suggests that turbulence plays an active role in the filament region. 

Finally comparison of both regions suggests a magnetic field interacting differently, and 
maybe less with its surrounding medium in the filament than it does in the core region, in 
agreement with the evolutionary stage of these two regions. Additional IR data would be 
necessary to confirm this point, and the quality of FIR measurements should of course be 
confirmed in both regions. 

4-2.2. From the filamentary clouds scale to the local Galactic scale 

The mean orientation of the ambient magnetic field can also be compared to the local 
Galactic magnetic field orientation. Figure 8 shows a polarization map in Galactic coordi- 
nates based on measurements compiled by Heiles (2000) with data for which Pjop > 3. The 
region presents a complicated picture because we are looking down a spiral arm and possibly 
because of the influence of the Cygnus star- forming complex. In their works Heiles (1996) 
and Walawender et al. (2001) have used a large set of polarization data and studied the 
inferred magnetic field geometry in this part of the Galaxy. At first sight, the simple picture 
of a toroidal magnetic field implies that the magnetic field lines should point in the direction 
toward / = 90°, b = 0°. In their statistical study Walawender et al. (2001) show that this 
trend is effectively real and are also able to detect the curvature of the local spiral arm. 
As mentionned by Jones (2003), since GF 9 is at a relatively close distance it is probably 
threaded by the magnetic field in the local spiral arm. Jones compares the orientation of 
the mean IR vector in the core region with the visible vectors on a greater scale, and shows 
that this vector points in the direction of the 'vanishing point' approximately located at 
I = 90°, b = 0°. We can see this vector denoted by 'IR' in Figure 8. We also show the mean 
visible vector 'V in this map. It does not point in the same direction than the 'IR' vector 
but it is approximately parallel to its nearest neighbours. All these facts suggest that it is 
preferable to use data probing the diffuse parts of the ISM in the environments of molecular 
clouds to compare ambient to the cloud and Galactic magnetic field's orientations. 
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4.3. Magnetic field strength and magnetic flux in the core region 

In order to estimate if the core region observed in CS by Ciardi et al. (2000) is 
magnetically supercritical or magnetically subcritical, we use the Chandrasekhar & Fermi 
(1953) method to estimate the magnetic field strength. The magnetic field strength in the 
plane-of-sky is estimated with the expression: 

/ A(V) 

B pos « 9.3v/^(Hj)^^G, (2) 

where 5<p is the polarization P.A. dispersion in degrees, n(H2) is the molecular hydrogen 
density in molecules cm~ 3 , AV = v8hi25(V) is the FWHM line width in km s _1 , and 5V is 
the velocity dispersion (see Crutcher (2004)). 

We use only the IR data for the core region since we have no way of estimating the 
turbulent velocity dispersion for the larger area corresponding to the visible data. The value 
5(f) = 5.9° is taken from Table 4 using the Jones IR data. This value is provided by a relatively 
small set of six measurements for which the mean of the uncertainties is 6.7°. If instead of 
a direct average, we use a weighted average for estimating the dispersion in P.A., we get 
7.2°. Finally, working with the P.A. directly instead of the vectors, the dispersions are 6.4° 
and 6.0°, weighted and unweighted respectively. All these values are close to each other and 
reflect the fact that there is little dispersion in P.A. in the cloud. Therefore, we believe that 
despite the small number of measurements, we nevertheless have a representative value for 
the dispersion. The various parameters that we used for the CS core region and our results 
are shown in Table 6. 

Once -Bpos is derived, we estimate A ODS , the observed mass to flux ratio (j^) observed 
normalized to its critical value. This observed normalized mass to flux ratio is subject to 
geometrical bias, as explained by Crutcher (2004). In the case of estimating the total 
field from the plane-of-sky component of the field, two of the three components are known 
and the correction factor for A Q bs is then 7r/4, which corresponds to A c given in Table 6. 
This correction is statistical in nature since we do not know the appropriate angle for any 
individual object. In fact, there may not be a line of sight component to the field. Using 5<p = 
8.9° for the dispersion in P.A. 3 to take into account the uncertainties in the measurements 
would decrease the magnitude of -B pos and increase the magnitude of A ODS and A c by about 
50%. For comparison, the dispersion in P.A. for the visible data around the core region 
is 12.4° (Table 4). With this value (even though we are outside the area for which the 
CS line-widths have been measured), the values of A Q bs and A c are larger by about 110%. 



3 8.9° is the quadratic sum of the dispersion 5.9° and the mean uncertainty 6.7° 
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These values are probably too large because the visible data cover a much larger area than 
just the core region. Also, a bias exists in the Chandrasekhar-Fermi method that tends to 
cause estimated field strentghs to be too large and thus to underestimate values of A Q b s (e.g. 
Heitsch et al. (2001)). 

What is the effect of the distance d on these results? The column density N is used to 
get the gas mass, so M oc r 2 , where r is the radius of the cloud, assumed to be spherical. 
The density n oc M / volume oc r" 1 . Since r oc d^ 1 , n oc d, B pos oc d 1 ^ 2 and A oc d^ 1 ! 2 . If the 
distance is 200 pc instead of 440 pc as assumed for the values given in Table 6, A would be 
closer to unity by about 

Is there additional possible evidence for a subcritical core? Crutcher (2004) identified 
two other tests to distinguish between ambipolar diffusion and turbulence as supporting 
mechanisms for star formation. If ambipolar diffusion is acting, one would expect an hour- 
glass configuration for the magnetic field. But the distribution of the IR polarization vectors 
(see Fig 6) suggests a smooth field, or even a slight curvature in the opposite direction than 
expected for an hourglass shape. A third test identified by Crutcher (2004), if ambipolar 
diffusion is acting one expects B oc p K with k approaching 1/2 for larger densitites (> 100 
cm -3 ). The values derived here are more or less compatible with this relation (see Figure 1 
in Crutcher (1999)). In addition, the comparison of the magnetic field orientations with the 
elongations of the cores seen in Table 5 is not in general compatible with ambipolar diffusion, 
since in that case the cores would be perpendicular to the direction of the magnetic field. 

Taking all the above arguments into account, our data suggest that the GF 9 core 
region is approximately magnetically critical. Ambipolar diffusion does not appear to be 
supported. Also, turbulence would produce a randomness in the polarization vectors which 
is not observed. The picture presented in section 4.2.1 of an original poloidal field which is 
transformed into a toroidal field by rotation appears most plausible. 

5. Conclusions 

Visible linear polarization data taken at the Mont-Megantic Observatory on 78 stars in 
the vicinity of the dark filamentary cloud GF 9 were presented and compared with available 
NIR and FIR polarization data covering part of the same region. Our main conclusions are: 

On a scale of a few times the mean radial dimension of the molecular clouds, the magnetic 
field is smooth and well-ordered, and appears to be 'blind' to the spatial distribution of the 
filaments. 
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On smaller scales in the core region, NIR polarimetry shows a rotation of the magnetic 
field lines when reaching into the denser parts. 

There is no evidence for saturation of polarization as a function of extinction for both the 
visible and NIR data, suggesting that even in the denser regions the same grain properties 
and grain alignment mechanisms apply, and that the magnetic field is still being probed 
adequately there. 

Finally, the magnetic field strength was estimated with the Chandrasekhar & Fermi 
(1953) method. The results would suggest that the core region is approximately magnetically 
critical. 

All these facts put together suggest that in GF 9 the magnetic field changes direction 
with the spatial scale. The field differs by 35° in the GF 9 core compared to the surrounding 
region. Ambipolar diffusion does not appear to be supported given the current evidence. 
The effects expected from turbulence are definitely missing. A global interpretation of the 
results is that an original poloidal field could have been twisted into a toroidal configuration 
by a rotating elongated structure during its collapse. 

Future work should check if the trend in polarization position angles to rotate coun- 
terclockwise with increasing wavelength still continues in the submm. Also, as the lack of 
saturation of polarization as a function of extinction suggests, if dust grains are magnetically 
aligned there should be no polarization hole in the submm in the core and the filament re- 
gions. As an additional benefit, submm polarization would afford a higher spatial resolution 
than the FIR data discussed here. 

We thank the Conseil de recherche en sciences naturelles et en genie du Canada for 
supporting this research. The authors thank J. S. Greaves for providing information about 
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We also thank an anonymous referee for his constructive comments and helpful suggestions. 
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Table 1. Visible polarization data. 
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see Table 2 
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Table 1 — Continued 
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1 


no 


65 


20:54:28.8 


+60:26:53.5 


1.42+0.21 


143+4.2 


2.32 


1 


0424700062 


62 


20:54:31.3 


+60:32:40.7 


2.53+0.15 


141.3+1.7 


1.90 


1 


no 


64 


20:54:53.0 


+60:29:35.7 


0.30+0.18 


102.7+17.2 


2.28 





0424700686 


73 


20:56:15.6 


+60:17:41.9 


1.89+0.15 


139.3+2.2 


1.94 


1 



( a 'l is tabulated when P > 3<rp, otherwise is tabulated. 
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Table 2. Stars with suspected variable polarization. Coordinates are given in Table 1. 



GF9 


P±a P 




J.D. 


number 


(%) 


(°) 


(days) 


35 


2.75±0.12 


135.8±1.3 


2,452,054.64583 


35 


0.43±0.13 


163.2±8.7 


2,452,141.72222 


86 


10.70±0.41 


131.3±1.1 


2,452,830.75347 


86 


12.41±0.37 


15.4±0.8 


2,452,835.78125 
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Table 3. NIR polarization data of Jones (2003) and equivalent A = 7660 A polarization 
data in the core and in the filament regions. Visible extinction coefficients are from 

Dobashi et al. (2005). 



Star 


P H ± o- Ph 


-P766O ± cp 7660 


A v 




(%) 


(%) 


(Mag.) 


CI 


1.9±0.6 


3.7 ±1.2 


2.34 


C2 


1.4±0.3 


2.8 ±0.8 


2.34 


C3 


4.2±1.0 


8.3 ±2.0 


2.60 


C4 


1.7±0.3 


3.4 ±0.8 


2.34 


C5 


4.0±0.9 


7.9 ± 1.8 


2.60 


C6 


2.7±0.6 


5.3 ± 1.2 


2.34 


F2 


1.9±0.4 


3.7 ±0.8 


2.19 
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Table 4. Means and dispersions of polarization and position angles at various wavelengths 
and various scales (see figures 6, 7, 2 and 8 respectively). 



Core region Filament region GF 9 scale Galactic scale 

P 7 660,V ± Sp 7660 ,v ( % ) 2 - 17 ± °- 93 2 - 62 ± L27 2 - 22 ± 133 L44 ± 126 <6) 

76 6O,v ±Se 7660 , v (°) (c) 136.0 ±12.4 126.1 ±13.1 127.4 ±25.5 72.0 ± 43.3 

P H ±S , p H (%)' 2.65 ±1.10 1.9 ± 0.4 M 

e H ± S Bh (°) (c) 171.2 ± 5.9 137.0 ± 6.0 < d > 
-Pfir ± Sp^iR (%) 

<?fir + 90° ± 5e FIR (°)< c ) » 19 random 



< a )Our data at A = 7660A. 

( ft )30° X 30° map centered at position (l c = 95°, b c = +10°) (see Figure 8), but in the V 
band. 

( c ) Position angles are given in the equatorial frame measured east from north. 

( d ) Only one measurement is used at this wavelength (see Figure 7). 
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Table 5. Position angles of cores in GF 9 and offsets with the mean magnetic field 
position angles as measured in the three different passbands. 



Core q(2000) 


5(2000) 


P.A>) 


|A(9 766 ol 


|A0 H | 


|A6> F ir| 


Location' 6 ' 


name ( h mn s ) 


(° ' ") 


(°) 


(°) 


(°) 


(°) 




LM349< C > 20 49 6.9 


60 11 48 


47 


ps 72 


ps 90 


random 


F 


LM350< C ) 20 49 46.5 


60 15 40 


49 


» 70 


ps 88 


random 


F 


LM351< C ) 20 51 28.0 


60 18 33 


43 


ps 91 


ps 51 


ps 24 


C 


LM352< C ) 20 53 49.9 


60 11 27 












L1082C( d > 20 51 27.4 


60 19 00 


7 


pa 53 


ps 15 


ps 12 


C 


L1082A( d ) 20 53 29.6 


60 14 41 












L1082B( d ) 20 53 50.3 


60 09 47 


48 


? 


? 


? 





( a ) East of north. 

( b ) F: filament region, C: core region. 

( c ) Lee & Myers (1999). 
WBenson & Myers (1989). 
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Table 6. Magnetic field strength in CS cores. 



Parameter 


GF9 Core 

(a) 


n (H 2 cm- 3 )! 4 ' 


15000 ± 3000 


TV (H 2 cm- 2 )!') 9.9 10 21 


R(pc)W 


0.16 


M (Mq)^) 


15 ±3 


T K (K)«>) 


7.5 ±0.5 


AV (fem/s)( fe ) 


0.88± 0.20 


«0 (°) 


5.9 




170±56 


A obs (c) 


ks 0.44 


Ac 


0.35 



( a )Only IR data are used. 
WCiardi et al. (2000). 

< c >A obs = 7.6 x HT 21 ^, 
Crutcher (2004). 
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Fig. 1. — DSS red plate of GF 9 and its environment with reverse colors. Various regions 
appearing in catalogs are identified in the Figure. GLMP 1012 is a post AGB star (see 
text for explanations). Boxes are regions observed by Ciardi et al. (1998), Ciardi et al. 
(2000) with letters C and F indicating the core region and the filament regions respectively. 
Positions of stars SI, S2, S15, S35, S70, S75, S86 and S92 are shown. The 'constellation' 
formed by stars S75, S70, S86, S35 and S92 can be used to compare the polarization pattern 
shown in Figure 2 with the location of the densest regions of the ISM as seen in the visible. 
The great black circle suggests the circular shape of the GF 9 filaments. The densest parts 
of the ISM forming GF 9 are located to the north of this circle, and small faint clouds are 
located in the south. 
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Fig. 2. — Visible polarization map of GF 9. The C and F rectangles are the core and the 
filament regions observed by Ciardi et al. (1998), Ciardi et al. (2000). Stars S75 and S92 
are embedded in the filament and have a polarization with a different orientation from the 
mean orientation of the ambient magnetic field. Object S70 is a visual binary. Crosses show 
the positions of stars S35 and S86 for which variable polarization is suspected (See Table 2). 
Only data for which P > 3o> are represented in the Figure. Bold vectors with sizes of 2% 
show position angles of cores LM349, LM350, LM 351, L1082C and L1082B. (See Figure 1). 
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Fig. 3. — Histogram of polarization in GF 9. Data for which P > 3a P (See last column in 
Table 1) are shown with bold lines. For this set of data, the mean and the dispersion of the 
distribution of the degrees of polarization are P = 2.22% and Sp = 1.33% respectively. 
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Fig. 4. — Histogram of position angles of polarization in GF 9. Data for which P > 3<7p 
(meaning oq < 9.5°) are shown with bold lines. For this set of data, the mean and the dis- 
persion of the distribution of the position angles are 9 = 127.4° and S$ = 25.5° respectively. 
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Fig. 5. — Variations of P with Eb-v- Ay data from Dobashi et al. (2005) were used to 
compute -Eb-v- The full line delineates the upper envelope P < 9.0P B -v (Serkowski et 
al. (1975)). A linear fit to the whole data set is shown by the dashed line and is given 
by P m 4.3-Eb-v- Dotted lines reproduce the truncation in the (P, E B -v) relation observed 
by Arce et al. (1998) in dark clouds in Taurus. X symbols: broadband red data from 
Observatoire du Mont-Megantic. Diamonds: NIR data from Jones (2003) converted to 
red-equivalent data (see section 3.3). 



-31 - 




Y09 £'09 Z'09 L '09 09 

(0002)-39Q 



Fig. 6. — Polarization map around and into the core (C) region observed by Ciardi et al. 
(1998), Ciardi et al. (2000). Bold vectors represent visible data. Thin vectors are IR data 
from Jones (2003) with P = 1.5% for clarity. Vectors with intermediate size are ISO FIR-90 
rotated data from Clemens et al. (1999) with P = 2.5% for clarity. Only FIR data for which 
the position angle is approximately the same when using the KLC intrumental polarization 
and the C-off instrumental polarization are shown. 
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Fig. 7. — Polarization map in the filament region. Same as in Figure 6, but around and into 
the filament (F) region observed by Ciardi et al. (1998), Ciardi et al. (2000). 
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Fig. 8. — Galactic polarization map around GF 9. Visible data from the Heiles catalog in a 
30° x 30° box centered at position (l c = 95°, b c = +10°) are shown. The IR symbol denotes 
the mean IR polarization of data taken in the GF9 core region by Jones (2003). The V 
symbol shows the mean visible polarization from our data in the vicinity of the core region. 



